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Abstract 

It is shown that particles undergoing discrete-time jumps in 3D, starting at a distance ro from the center of 
an adsorbing sphere of radius R, are captured with probability (R — ca)/ro for ro 3> R, where c is related 
to the Fourier transform of the scaled jump distribution and a is the distribution's root-mean square jump 
length. For particles starting on the surface of the sphere, the asymptotic survival probability is non-zero 
(in contrast to the case of Brownian diffusion) and has a universal behavior a/ (Ry/E) depending only upon 
a/R. These results have applications to computer simulations of reaction and aggregation. 
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I. INTRODUCTION 



A celebrated result in probability theory states that in three or higher dimensional space, par- 
ticles undergoing diffusion or Brownian motion have a chance to escape to infinity rather than be 
captured to a finite adsorbing surface. This was first demonstrated in the case of random walks 
(discrete versions of Brownian motion) on a 3D cubic lattice, on which a walker starting at a site 
adjacent to the adsorbing origin escapes with probability J -1 « 0.659463 or is captured with 
probability 1 — J -1 0.340537, where flfl ^] 

dx dy dz 



(27r) 3 J_ n J_ w J_ w 3 — cosx — cosy — cos z 

= ^r(£)r(£)r(£)r(£) . (i) 

In fewer than three dimensions, walks return to the origin an infinite number of times and are thus 
recurrent; while in three and higher dimensions they escape and are transient. The transient nature 
of random walks in 3D is the basis for long-range diffusive transport and is essential for many 
important physical phenomena involving particles that must escape a locally trapping region. 

The same phenomenon in continuum theory is illustrated by the problem of adsorption of dif- 
fusing point particles on a spherical boundary of radius R. Say the particle starts initially at 
a distance |r | > R from the center of a sphere. The survival probability S(r ,t) up to time 
t satisfies the backward diffusion equation dS(r ,t)/dt = DV 2 S(ro,t) with the initial con- 
dition S(r ,0) = 1 for all |r | > R and the boundary conditions S(\r \ = R,t) = and 
S(\r \ — > oo,t) = 1 for all t. Exploiting radial symmetry, the solution to this equation is eas- 
ily found at all times t 

S(r ,t) = l--erfc r -^£. (2) 
Taking the limit t — > oo, one finds that the probability to ultimately escape is simply 

S(r , oo) = 1 - R/r , (3) 

a well-known result in probability theory. The capture probability is P(r , oo) = 1 — S(r , oo) = 
R/r . 

In this paper, we ask how having a finite step-length in the diffusing particle's motion and dis- 
crete time-steps affects the continuum results in Eqs. (1213b . In particular, we consider a random 
walker in 3D continuum where, in each step of equal time interval r, the particle jumps isotrop- 
ically a distance r drawn from a distribution A7rr 2 W(r), bounded above by 2R, but otherwise 



FIG. 1 : Example of a walk that starts near the surface of the sphere of radius R and eventually escapes. The 
dashed circle represents the effective capture radius R' = R — ca for particles that start far from the sphere. 



arbitrary. If the particle jumps inside the sphere it is captured. Let S n (r ) be the probability 
that the particle, starting at r > R, survives up to step n. The net capture probability is then 
-Poo( r o) — 1 — 5oo( r o)- The question of the net capture probability in this case is also important to 
many computer simulation problems, where particle diffusion is carried out by such a discrete-time 
process, typically through a Pearson flight in which the particle at each time step jumps a discrete 
distance £ which may be of the order of magnitude of R or even larger. The capture probability can 
be used to find chemical or catalytic reaction rates, and is also relevant to non-equilibrium growth 
studies, notably in diffusion-limited aggregation [4] and related aggregation models. There has 
been little (or no) discussion on how the simulation's jump length affects the ultimate results of 
these simulations. 

Another situation where a finite jump length is relevant is when particles are diffusing in a gas 
and the mean free path becomes comparable to the diameter of the adsorbing sphere. Then the 
motion of the diffusing particle is not an infinitesimal Wiener process described by the diffusion 
equation, but instead is described by the Ornstein-Uhlenbeck theory of Brownian motion with 
inertial drift, integrated to position space [|5|]. Only after particles travel a certain distance do 
they "forget" their initial direction and continue in a random direction. An approximation to this 
process is given by discrete random walk. 

For the discrete random walk problem, we find that the survival probability S n (r ) becomes 
substantially modified from its continuous-time counterpart S(r ,t). We show that the ultimate 
survival probability Soc(r ) has the following asymptotic behavior 

SooOo) ~ 1 for r > R 

~ m as ro " R (4) 

where a 2 = An J W(r)r 4 dr is the mean-square jump distribution and c is a constant (computed 
explicitly below) that depends upon the details of W(r). Thus if the particle starts very far from the 



surface of the sphere, the final survival probability has a similar expression as in the continuous- 
time case in Eq. ©, except that the effective radius of the sphere R' = R — co shrinks by a 
finite 'extrapolation' length ca (see Fig. [T]). On the other hand, if the particle starts right outside 
the surface of the sphere we find a rather surprising result: in sharp contrast to the continuous- 
time case, the survival probability is nonzero and moreover has a universal value 1 / \/6 in units 
of the dimensionless ratio a/R, and is completely independent of the other details of the jump 
distribution (see Fig. |2]). 



II. DERIVATION OF RESULTS 



To derive these results our strategy is to first map the capture problem to the so called 'flux' 
problem to a sphere and then use the known results of the latter problem. This mapping works 
both for continuous-time Brownian motions as well as for discrete-time random walkers. For 
simplicity we discuss first the mapping for continuous-time Brownian motions. The continuous- 
time flux problem, first studied by Smoluchowski [(], 7], is defined as follows. Consider an infinite 
number of noninteracting Brownian particles initially distributed with uniform density p outside a 
sphere of radius R in 3D. Each particle subsequently performs Brownian motion and, if it reaches 
the sphere, it is absorbed. If p(r, t) denotes the density profile at time t, the instantaneous flux to 
the sphere is = 4iiR 2 Dd r p(r, t)\ r= R. The density profile p(r, t) can be calculated easily as 
follows. For the uniform initial condition, p(r, t) is only a function of the radial distance r and 
it satisfies the diffusion equation d t p = D [d^p + (2/r)<9 r p] for r > R with the initial condition 
p(r, 0) = p an d the boundary conditions p(R,t) = and p(oo,t) = p . One can reduce this 
problem to a ID diffusion problem with F(r, t) = rp(r, t) satisfying d t F(r, t) = Dd^F(r, t), and 
F(r, 0) = p r and F(r = R,t) = 0, whose solution can be easily found via the method of images. 
Dividing F(r, t) by r one gets the density profile 



p(r, t) = Po 



R r r-R 

— erfc— - r 

r ^/Wi 



(5) 



Consequently the flux is = AuRDpoyl + R/y/nDt] which tends to a constant $(oo) = 
AixRDpQ at long times. Comparing with Eq. © one sees that the survival probability S(r , t) has 
the same expression as the density profile p(r, t)/po in the flux problem provided one replaces r 
by r in the latter problem. 

To understand the origin of this connection between the two problems it is useful to discuss both 
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of them simultaneously in terms of a single Green's function. The Green's function G(r, r , t) for 
finding the particle at position r at time t, starting at r outside the sphere at t = 0, satisfies 

dG(r,T ,t) nv7 2^/ .X ,~ 
— = DV G(r, r , t) , (6) 

subject to the initial condition G(r,r o ,0) = 6(r — r ), and the absorbing boundary condition 
G(t, r , t) — for |r| = R. In the flux problem, given an arbitrary initial density p(r , 0), the 
density at time t can be found from 

p(r,t) = J G(r,r ,t)p(r ,0)dr (7) 

where the integration is over r > R. In the capture problem, the survival probability S(tq, t) up 
to time t, starting at r , can also be written in terms of the same Green's function, 



S(r ,t) = J G(r,r ,t)dr 



(8) 



where the integration is over r > R. For uniform and isotropic initial density p m the flux 
problem one gets from Eq. © 



p(r,t) = p y G(r,r ,t)dr . 



(9) 



Since the Green's function G(t, r , t) is symmetric in r and r for unbiased Brownian motion, it 
follows by comparing Eqs. ([8]) and © the equality 

S(r Q ,t)=p(T ,t)/po (10) 

This leads to the quite general result that the probability 5(r , t) a particle at position r survives 
up to time t can be found directly by solving the density profile in the corresponding flux problem 
starting from an initial uniform density. Eq. (flOl) just says that the net probability a walker starting 
at point r survives up to time t is the same as the sum of the probabilities that particles starting 
at every point in space reach r in time t without being adsorbed. We note that the same general 
argument goes through even for discrete-time random walks, leading to the relation S n (ro) = 
p n (r )/p , where p„(r ) is the density profile at step n of the corresponding discrete-time flux 
problem. This equivalence is also discussed in Refs. [8, 
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Thus we can use many of the discrete-time flux-problem results from Refs. ^[kIiUI] to study 



the present capture problem. Besides mapping those results onto the present problem, we derive 



several new results, details of which will be given in II 1211 . We also reformulate the mathematical 



expressions in order to scale out the root-mean square (r.m.s.) jump length o from the integral 



expressions and leave them dimensionless ft 1 311 - 

We consider that a particle at position r' (with r' > R) jumps in one time step (of duration r) to 
a new position r with a jump distance |r — r'| that is drawn independently from an isotropic distri- 
bution 4irr 2 W(\r — r'|), bounded above by 2R, and normalized as Air J^ R W(r)r 2 dr = 1, as illus- 
trated in Fig.Q] Because of radial symmetry, this problem can be formulated as a one-dimensional 
problem, exactly analogous to the case of radial diffusion discussed above, with effective density 
p n (r) = rp n (r) and an effective jump probability given by the symmetric, non-negative function 

f(x) = 2vr / W(r)rdr . (11) 

J\x\ 

which is also normalized to unity, f(x)dx = 1. The r.m.s. jump length of / is l/y/3 times the 
jump length for W: a 2 = J °° W(r) r 2 Anr 2 dr = 3 f(x) x 2 dx. The quantity p n (r) satisfies 

rr+lm 

Pn+i( r ) = / f(\r -r'\)p n (r')dr' (12) 

J\r-i m \ 

where £ m is the maximum of the jump distance, and the subscript n represents the time step. 

We rescale the jump distribution by a and thus define a new function g(y) by f(x) = 
(l/a)g(x/a) so that g (y) is normalized to unity and has a second moment of 1/3. Two important 
quantities related to g(y) appear: 
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oo 

In 



i-g(k) 



and 



ttV6 Jo 

where g(k) = g(y)e lky dy. 



k 2 /Q 

POO 

/ m[l — dk (14) 
Jo 



In llllh . the general solution to Eq. (fT2~l) is given, in terms of a double Laplace transform of 
S n (f) — Pn( r ) I Po- Here we analyze that result explicitly in two important limiting cases: 

(i) For r 3> R and for a large number of time-steps n 3> 1, the discrete-time survival proba- 
bility S n (r ) behaves as 

Sn(ro) = 1 - ^erfc^=4= + 0{n^ 2 ) (15) 

where R' = R — ca. This is valid for r <^ na, in which case r is well within the maximum 
distance a particle can travel in n bounded steps. We can write the above result explicitly in terms 



of time t = nr, and introduce the effective diffusion coefficient 

D = a 2 /(Qr) (16) 

which implies that the factor 2a 2 n/3 in (TT3T ) equals ADt. Then Eq. (fT3T ) becomes identical to 
Eq. (O, except that the effective radius of the capture sphere is reduced, as shown in Fig. [2] The 
ultimate survival probability for r 3> R is simply Soo(r ) = 1 — R' /r as given in Eq. © 

Note that if we take the Brownian limit of a — > and r — >■ with D = const., then Eq. (fl~5l) 
becomes Eq. © exactly. 

(ii) For n = oo, the Laplace transform of the steady-state solution, written in terms of g(k), 
simplifies to 

F^e-^dz 

A /"» ln[l - $(*)] 



\V6 



exp 



-dk 



(17) 



Wo A 2 + fc 2 

where F 00 (z) = r S (XJ (r )/a and z = (r — i?)/cr. For z <C 1, Eq. (TTTI) implies ^00(2) = 
1 / v^6 + . . ., which yields 

&oW = ' + (18) 

r V6 ^0 

for r — R cr. When r = -R, this gives Eq. ©. If we set a = R (the mean jump length equal 
to the radius), the escape probability from the surface is 1/ y/E ~ 0.408248. This compares with 
the escape probability w 0.659463 for a cubic lattice that follows from Eq. (OQ). 

It is interesting to note that as a — >• 0, the first term in Eq. (fl~8~l) vanishes, but the second term 
does not. Thus, no matter how small a/R is, in the small region near the surface R < r <^ R + a, 
the slope of the curve of (r ) vs. r has the value rather than the value 1/R corresponding 
to the diffusion-equation solution. 

Because the asymptotic flux $(£—>• 00) has a universal value AixDRp Q in the Brownian limit, 
one might wonder whether the universality of (R) given by Eq. © (now viewed in the context 
of the flux problem) is related to that universal flux value. The flux $ n as n — ► 00 is found from 



L12] 



R+e m rR 



$^ = 1^/ dr'r'S^r') I drrf(r-r') (19) 

T JR Jr'-i m 

and depends upon the behavior of S^ir) within a distance i m of the sphere. For to give 
A%DRp to leading order in a/R, it follows that we must have 

dz'F^z') f dzg(z-z') = I (20) 
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FIG. 2: Behavior of the net (infinite-time) survival probability S^ro) as a function of the distance r$ from 
the center of the adsorbing sphere (solid line). Also shown is the continuous-time prediction (heavy dashed 
line), and extrapolation of asymptotic behavior to = (light dashed line). The extrapolation length is 
ca, the actual value of Soo(R) is (r/(Ry6), while its extrapolated value is ca/R. 

But all orders in the expansion of F oa (z') = l/v6 + bz' ... contribute to the integral in Eq. d20l) . 
and so it follows that the universal leading term in Eq. (fl~8l) is not simply a consequence of the 
universality of the flux, making that term all the more intriguing. 

The values of c and b can be calculated for various jump distributions. The most natural one 
from the point of view of computer simulation is the Pearson flight W(r) = 6(r — £)/(4:Ti£ 2 ), 
which leads to f(x) = 1/(2/) for x £ (—£, I) or g(y) = 1/2 for y £ (—1, 1) (and zero otherwise), 
andg(A;) = smk/k. Then 



c=-i f.n 

Jo 



(sinO/e 



£76 



s c2 



-± « 0.2979521903 



and 



7T 



v/6 



In 



sin£ 



d£ « 0.6538250956 



(21) 



(22) 



The number c = 0.29795 . . . has had a fairly long history. It first appeared in [|9|] in the context 
of the flux problem, where its value was determined through a numerical iteration. It later appeared 
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independently 111411 in a study of the asymptotics of a sum of random variables with a uniform 
distribution, and was evaluated by a slowly converging double summation. The integral form (l2"TT) 
was given in [15] in the context of the ID random-walk problem. Finally, relation between the 3D 
flux and ID random-walk problems was shown in I loll . 

For the Pearson flight, we can find F^z) to higher order numerically F^z) = l/y/E + 
bz + 0.2362658938^ 2 + 0.014221827913,^ . . ., and in this case we confirm that Eq. (f20]) is indeed 
satisfied to high precision (when more of these terms are included). For the Pearson flight, the 
second integral in Eq. (|20l) gives (1 — z')/2. Also, if we define Ik = f (1 — z) h F 00 (z)dz, then we 
find I = 2/y/E, h = 1/3, and J 2 = 2c/3. We also find 2^,(1) = 2b. 

As a second example, we consider that the jump is uniform within a sphere of radius I, so that 
W(r) = 3/(4n£ 3 ) for |r| < I and zero otherwise, implying f(x) = (3/4)(£ 2 — x 2 )/£ 3 , g(y) = 
(3/4)v^75(l-3y 2 /5)(|y| < y/5/3Xg(k) = 9VTEsm(k^/5/3)/(25k 3 ) -9cos(k^/5/3)/(5k 2 ), 
b = 0.682012 ... and c = 0.310901 . . .. Notice that the value of c is just a little larger than that of 
the Pearson walk. 



III. FURTHER DISCUSSION AND CONCLUSIONS 



In deriving these results, we assumed that adsorption occurs only if the final position of the 
particle falls within the sphere - thus trajectories that pass through or graze the sphere are assumed 
not to adsorb. Including these events lead to corrections of higher order in a/R flSfl. Physically, 
of course, such trajectories should be adsorbed, so in that case our results are only accurate for 
a/R 1 . However, for computer simulation it is easiest and therefore common to only check for 
adsorption based upon the final position, in which case our results are exact as long as W(r) = 
for r > 2R. 

The assumption W(r) = for r > 2R is necessary so that the transformation from the 3D 
to the ID problem described by (|T2|) is exact. Once on the ID level, however, the mathematics 
of the calculation holds for any length jump distribution, and here we consider two infinite-range 
models. If a R, then the probability of jumping beyond 2R is very low and these results should 
give good approximations to the true 3D behavior. 

In the first of these models, we consider the exponential distribution g(y) = a/3/2 e^^ y \ 
implying^) = 6/(k 2 + 6), and find c = 1 / V® = 0.408248... and b = 1. Thus, R' = R-c/y/E, 
and by Eq. (fT8l) , S^Vq) = a/(r Q y/6) + (r — R)/r = (r — R + ca)/r , which in this case is the 
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solution for all r , not just for r — R <C o. 

Secondly, we consider a distribution that is Gaussian in W and therefore also in g. Then, 
g(y) = y/3/(2n) exp(-3y 2 /2\_g(k) = exp(-k 2 /6), and c = -t(l/2)/y/6n = 0.336363... 
(similar to what was found in jl5\\ in another context) and b = £(3/2)/(2-y/^r) = 0.736937. . .. 
Interestingly, C(l/2) a l so appears in the problem of adsorption of a particle (in ID) diffusing by 
the Ornstein-Uhlenbeck process . 11811 . where it is also related to the boundary extrapolation 



length, as well as in in the context of the maximum of a Rouse polymer chain lil5h or Gaussian 
random walks J19I1 . 

Thus, we have shown that the discreteness of the time always affects the ultimate capture prob- 
ability of a particle undergoing a random walk. Even if the walk length were drawn from a proper 
Gaussian distribution appropriate for Brownian motion representing the time interval r, the cap- 
ture probability would not be the same as the solution for the diffusion equation, because of the 
ability of the discrete walk to jump away from the surface. Different distributions change the 
constants c and b, but the general behavior remains the same. 

This work implies that, when doing computer simulations involving capture on a sphere, in or- 
der to get a capture probability that is say 99% correct, one should make a less than (0.01/c)i2, or 
in the case of the Pearson flight, make i less than 0.03356i? by Eq. (f2Tb . Indeed, most simulations 
use much larger jump lengths and therefore significantly underestimate the capture probability. 
An improvement in accuracy can be achieved by varying the step size as the adsorbing boundary 



is approached ll2D. 
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